Validation of a forage production index (FPI) derived from MODIS fcover time-series using high-resolution satellite imagery: methodology, results and opportunities

An index-based insurance solution was developed to estimate and monitor near real-time forage production using the indicator Forage Production Index (FPI) as a surrogate of the grassland production. The FPI corresponds to the integral of the fraction of green vegetation cover derived from moderate spatial resolution time series images and was calculated at the 6 km x 6 km scale. An upscaled approach based on direct validation was used that compared FPI with field-collected biomass data and high spatial resolution (HR) time series images. The experimental site was located in the Lot and Aveyron departments of southwestern France. Data collected included biomass ground measurements from grassland plots at 28 farms for the years 2012, 2013 and 2014 and HR images covering the Lot department in 2013 (n = 26) and 2014 (n = 22). Direct comparison with ground-measured yield led to good accuracy (R-2 = 0.71 and RMSE = 14.5%). With indirect comparison, the relationship was still strong (R-2 ranging from 0.78 to 0.93) and informative. These results highlight the effect of disaggregation, the grassland sampling rate, and irregularity of image acquisition in the HR time series. In advance of Sentinel-2, this study provides valuable information on the strengths and weaknesses of a potential index-based insurance product from HR time series images

Potential improvement for forest cover and forest degradation mapping with the forthcoming sentinel-2 program

The forthcoming European Space Agency’s Sentinel-2 mission promises to provide high (10 m) resolution optical data at higher temporal frequencies (5 day revisit with two operational satellites) than previously available. CNES, the French national space agency, launched a program in 2013, ‘SPOT4 take 5’, to simulate such a dataflow using the SPOT HRV sensor, which has similar spectral characteristics to the Sentinel sensor, but lower (20m) spatial resolution. Such data flow enables the analysis of the satellite images using temporal analysis, an approach previously restricted to lower spatial resolution sensors. We acquired 23 such images over Tanzania for the period from February to June 2013. The data were analysed with aim of discriminating between different forest cover percentages for landscape units of 0.5 ha over a site characterised by deciduous intact and degraded forests. The SPOT data were processed by one extracting temporal vegetation indices. We assessed the impact of the high acquisition rate with respect to the current rate of one image every 16 days. Validation data, giving the percentage of forest canopy cover in each land unit were provided by very high resolution satellite data. Results show that using the full temporal series it is possible to discriminate between forest units with differences of more than 40% tree cover or more. Classification errors fell exclusively into the adjacent forest canopy cover class of 20% or less. The analyses show that forestation mapping and degradation monitoring will be substantially improved with the Sentinel-2 programJRC.H.3-Forest Resources and Climat

Irrigated grassland monitoring using a time series of terraSAR-X and COSMO-skyMed X-Band SAR Data

[Departement_IRSTEA]Territoires [TR1_IRSTEA]SYNERGIE [Axe_IRSTEA]TETIS-ATTOSInternational audienceThe objective of this study was to analyze the sensitivity of radar signals in the X-band in irrigated grassland conditions. The backscattered radar signals were analyzed according to soil moisture and vegetation parameters using linear regression models. A time series of radar (TerraSAR-X and COSMO-SkyMed) and optical (SPOT and LANDSAT) images was acquired at a high temporal frequency in 2013 over a small agricultural region in southeastern France. Ground measurements were conducted simultaneously with the satellite data acquisitions during several grassland growing cycles to monitor the evolution of the soil and vegetation characteristics. The comparison between the Normalized Difference Vegetation Index (NDVI) computed from optical images and the in situ Leaf Area Index (LAI) showed a logarithmic relationship with a greater scattering for the dates corresponding to vegetation well developed before the harvest. The correlation between the NDVI and the vegetation parameters (LAI, vegetation height, biomass, and vegetation water content) was high at the beginning of the growth cycle. This correlation became insensitive at a certain threshold corresponding to high vegetation (LAI ~2.5 m2/m2). Results showed that the radar signal depends on variations in soil moisture, with a higher sensitivity to soil moisture for biomass lower than 1 kg/m². HH and HV polarizations had approximately similar sensitivities to soil moisture. The penetration depth of the radar wave in the X-band was high, even for dense and high vegetation; flooded areas were visible in the images with higher detection potential in HH polarization than in HV polarization, even for vegetation heights reaching 1 m. Lower sensitivity was observed at the X-band between the radar signal and the vegetation parameters with very limited potential of the X-band to monitor grassland growth. These results showed that it is possible to track gravity irrigation and soil moisture variations from SAR X-band images acquired at high spatial resolution (an incidence angle near 30°)

A snow cover climatology for the Pyrenees from MODIS snow products

International audienceThe seasonal snow in the Pyrenees is critical for hydropower production, crop irrigation and tourism in France, Spain and Andorra. Complementary to in situ observations , satellite remote sensing is useful to monitor the effect of climate on the snow dynamics. The MODIS daily snow products (Terra/MOD10A1 and Aqua/MYD10A1) are widely used to generate snow cover climatologies, yet it is preferable to assess their accuracies prior to their use. Here, we use both in situ snow observations and remote sensing data to evaluate the MODIS snow products in the Pyrenees. First, we compare the MODIS products to in situ snow depth (SD) and snow water equivalent (SWE) measurements. We estimate the values of the SWE and SD best detection thresholds to 40 mm water equivalent (w.e.) and 150 mm, respectively , for both MOD10A1 and MYD10A1. κ coefficients are within 0.74 and 0.92 depending on the product and the variable for these thresholds. However, we also find a seasonal trend in the optimal SWE and SD thresholds, reflecting the hysteresis in the relationship between the depth of the snow-pack (or SWE) and its extent within a MODIS pixel. Then, a set of Landsat images is used to validate MOD10A1 and MYD10A1 for 157 dates between 2002 and 2010. The resulting accuracies are 97 % (κ = 0.85) for MOD10A1 and 96 % (κ = 0.81) for MYD10A1, which indicates a good agreement between both data sets. The effect of vegetation on the results is analyzed by filtering the forested areas using a land cover map. As expected, the accuracies decrease over the forests but the agreement remains acceptable (MOD10A1: 96 %, κ = 0.77; MYD10A1: 95 %, κ = 0.67). We conclude that MODIS snow products have a sufficient accuracy for hy-droclimate studies at the scale of the Pyrenees range. Using a gap-filling algorithm we generate a consistent snow cover climatology, which allows us to compute the mean monthly snow cover duration per elevation band and aspect classes. There is snow on the ground at least 50 % of the time above 1600 m between December and April. We finally analyze the snow patterns for the atypical winter 2011–2012. Snow cover duration anomalies reveal a deficient snowpack on the Span-ish side of the Pyrenees, which seems to have caused a drop in the national hydropower production

Kalideos OSR MiPy : un observatoire pour la recherche et la démonstration des applications de la télédétection à la gestion des territoires

International audienceCes dernières années, le CESBIO a mis en place un Observatoire Spatial Régional 'OSR' : un dispositif d'observation couplant mesures de terrain et télédétection dans le sud-ouest de la France. L'OSR se base sur des acquisitions mensuelles de données satellitaires à résolution décamétrique depuis 2002 et sur des sites expérimentaux lourdement instrumentés (mesures en continu de flux d'eau et de carbone) à partir de 2004. Ce dispositif a été reconnu service d'observation par l'INSU/CNRS en 2007 et site KALIDEOS par le CNES fin 2009 : 'KALIDEOS OSR MiPy'. Le site atelier correspond à une emprise d'image SPOT, soit environ 50x50 km et couvre une grande diversité de milieux (pédologie, topographie), d'occupation et d'utilisation des sols, de pratiques et de modalités de gestion (agricole, forestière...) et de conditions climatiques (fort gradient de déficits hydriques estivaux). Pour la télédétection, ce site a servi la préparation de SMOS, et il soutient maintenant en priorité à la préparation des missions VENμS et Sentinel-2. Les aspects radar, imagerie thermique et les approches multi-capteurs se développent depuis peu. Le traitement du signal, la physique de la mesure et l'amélioration de la qualité des données constituent le premier axe de recherche. Au niveau thématique, le CESBIO a pour priorité les suivis et les modélisations des agrosystèmes de grandes cultures. L'implication récente d'autres partenaires scientifiques ou gestionnaires a permis d'initier des travaux sur d'autres aspects, comme la biodiversité, l'aménagement du territoire, le suivi de l'extension urbaine, les risques environnementaux, la santé des forêts, l'enfrichement, la diversité et la productivité des prairies. La valorisation des 10 années d'archives 2002-2011 débute et semble très pertinente pour la caractérisation en haute et en basse résolution des conséquences d'années climatiques atypiques (2003, 2011) sur les éco-agro-systèmes. L'extrapolation des résultats obtenus sur ce site atelier à toute la région Midi-Pyrénées ou à la chaine des Pyrénées est aussi initiée

Uncertainty assessment of surface net radiation derived from Landsat images

The net radiation flux available at the Earth's surface drives evapotranspiration, photosynthesis and other physical and biological processes. The only cost-effective way to capture its spatial and temporal variability at regional and global scales is remote sensing. However, the accuracy of net radiation derived from remote sensing data has been evaluated up to now over a limited number of in situ measurements and ecosystems. This study aims at evaluating estimates and uncertainties on net radiation derived from Landsat-7 images depending on reliability of the input surface variables albedo, emissivity and surface temperature. The later includes the reliability of remote sensing information (spectral reflectances and top of canopy brightness temperature) and shortwave and longwave incoming radiations. Primary information describing the surface is derived from remote sensing observations. Surface albedo is estimated from spectral reflectances using a narrow-to-broadband conversion method. Land surface temperature is retrieved from top of canopy brightness temperature by accounting for land surface emissivity and reflection of atmospheric radiation; and emissivity is estimated using a relationship with a vegetation index and a spectral database of soil and plant canopy properties in the study area. The net radiation uncertainty is assessed using comparison with ground measurements over the Crau–Camargue and lower Rhone valley regions in France. We found Root Mean Square Errors between retrievals and field measurements of 0.25–0.33 (14–19%) for albedo, ~ 1.7 K for surface temperature and ~ 20 W·m− 2 (5%) for net radiation. Results show a substantial underestimation of Landsat-7 albedo (up to 0.024), particularly for estimates retrieved using the middle infrared, which could be due to different sources: the calibration of field sensors, the correction of radiometric signals from Landsat-7 or the differences in spectral bands with the sensors for which the models where originally derived, or the atmospheric corrections. We report a global uncertainty in net radiation of 40–100 W·m− 2 equally distributed over the shortwave and longwave radiation, which varies spatially and temporally depending on the land use and the time of year. In situ measurements of incoming shortwave and longwave radiation contribute the most to uncertainty in net radiation (10–40 W·m− 2 and 20–30 W·m− 2, respectively), followed by uncertainties in albedo (< 25 W·m− 2) and surface temperature (~ 8 W·m− 2). For the latter, the main factors were the uncertainties in top of canopy reflectances (< 10 W·m− 2) and brightness temperature (5–7 W·m− 2). The generalization of these results to other sensors and study regions could be considered, except for the emissivity if prior knowledge on its characterization is not available

The MODIS (collection V006) BRDF/albedo product MCD43D: temporal course evaluated over agricultural landscape

The assessment of uncertainties in satellite-derived global surface albedo products is a critical aspect for studying the climate, ecosystem change, hydrology or the Earth's radiant energy budget. However, it is challenged by the spatial scaling errors between satellite and field measurements. This study aims at evaluating the forthcoming MODerate Resolution Imaging Spectroradiometer (MODIS) (Collection V006) Bidirectional Reflectance Distribution Function (BRDF)/albedo product MCD43D over a Mediterranean agricultural area. Here, we present the results from the accuracy assessment of the MODIS blue-sky albedo. The analysis is based on collocated comparisons with higher spatial resolution estimates from Formosat-2 that were first evaluated against local in situ measurements. The inter-sensor comparison is achieved by taking into account the effective point spread function (PSF) for MODIS albedo, modeled as Gaussian functions in the North–South and East–West directions. The equivalent PSF is estimated by correlation analysis between MODIS albedo and Formosat-2 convolved albedo. Results show that it is 1.2 to 2.0 times larger in the East–West direction as compared to the North–South direction. We characterized the equivalent PSF by a full width at half maximum size of 1920 m in East–West, 1200 m in North–South. This provided a very good correlation between the products, showing absolute (relative) Root Mean Square Errors from 0.004 to 0.013 (2% to 7%), and almost no bias. By inspecting 1-km plots homogeneous in land cover type, we found poorer performances over rice and marshes (i.e., relative Root Mean Square Error of about 11% and 7%, and accuracy of 0.011 and − 0.008, respectively), and higher accuracy over dry and irrigated pastures, as well as orchards (i.e., relative uncertainty < 3.8% and accuracy < 0.003). The study demonstrates that neglecting the MODIS PSF when comparing the Formosat-2 albedo against the MODIS one induces an additional uncertainty up to 0.02 (10%) in albedo. The consistency between fine and coarse spatial resolution albedo estimates indicates the ability of the daily MCD43D product to reproduce reasonably well the dynamics of albedo

Agrometerological study of semi-arid areas : an experiment for analysing the potential of time series of FORMOSAT-2 images (Tensift-Marrakech plain)

Earth Observing Systems designed to provide both high spatial resolution (10m) and high capacity of time revisit (a few days) offer strong opportunities for the management of agricultural water resources. The FORMOSAT-2 satellite is the first and only satellite with the ability to provide daily high-resolution images over a particular area with constant viewing angles. As part of the SudMed project, one of the first time series of FORMOSAT-2 images has been acquired over the semi-arid Tensift-Marrakech plain. Along with these acquisitions, an experimental data set has been collected to monitor land-cover/land-use, soil characteristics, vegetation dynamics and surface fluxes. This paper presents a first analysis of the potential of these data for agrometerological study of semi-arid areas